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Using Z"f candidate events collected by the CDF detector at the Tevatron Collider, we search 
for potential anomalous (non-standard-model) couplings between the Z boson and the photon. Z'y 
couplings vanish at tree level and are heavily suppressed at higher orders; hence any evidence 
of couplings indicates new physics. Measurements are performed using data corresponding to an 
integrated luminosity of 4.9 fb~^ in the Z ^ vv decay channel and 5.1 fb~^ in the Z — >■ V^l~ 
{I — fi, e) decay channels. The combination of these measurements provides the most stringent 
limits to date on Z^ trilinear gauge couplings. Using an energy scale of A = 1.5 TeV to allow for 
a direct comparison with previous measurements, we find limits on the CP-conserving parameters 
that describe Zj couplings to be \hl"^ \ < 0.022 and | < 0.0009. These results are consistent 
with standard model predictions. 
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Studies of triliiicar couplings between the gauge bosons 
{W, Z, 7) test the standard model (SM) description 
of gauge sector interactions and provide sensitivity to 
physics beyond the SM through examination of produc- 
tion rates and kinematics [l|-l6|. In the case of neutral 
couplings, ZZ"f and Zjj vertex interactions vanish at 
tree level and, while allowed via internal particle loops, 
are highly suppressed in the SM. However, these trilinear 
gauge couplings can be non-negligible if loop contribu- 
tions occur via non-SM particles. Models such as those 
incorporating compositeness or supersymmetry can alter 
the predicted cross section and production kinematics of 
Z7 events 0-[l3l- 

In the SM, given the suppression of ZZj and ^77 cou- 
plings, the production of Z7 events is dominated by pro- 
duction of a Z boson along with the radiation of a pho- 
ton off either an incoming parton or a Z decay product. 
These production mechanisms are interesting in their own 
right, serving as an important background to searches 
for new physics (e^. in gauge-mediated supersymmetry 
breaking models 
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and Higgs boson searches. In this 
Letter, the production properties of Z7 events are com- 
pared to SM predictions, and limits are set on anomalous 
trilinear gauge couplings. 

The measurements of Zj couplings are performed with 
pp collision data at ^/s = 1.96 TeV from the Tevatron Col- 
lider using the Collider Detector at Fermilab (CDF). We 
seek two types of Z7 events: those where the Z decays 
to charged leptons (by identifying lepton candidate pairs 
and a p rompt photon [l^ with large transverse energy 
Ex [l3|), and those where the Z decays to neutrinos (by 



identifying an event with only a solitary, prompt, high- Et 
photon). In the former case, data corresponding to an in- 
tegrated luminosity of 5.1 fb"-'^ are used; in the latter, 4.9 
fb~^. These measurements use over twice as much data 
as the previous published CDF result [l[ and incorporate 
looser muon selection requirements. As no significant dis- 
agreement is found between the SM prediction and the 
data, we set limits that are not only far more restrictive 
than those measured in [l[, but are approximately half 
the magnitude of the previous best published limits Q ■ 

In beyond-the-SM scenarios with enhanced Zj cou- 
plings, not only does the Z7 production cross section 
increase, but the photon Et spectrum is modified due 
to an enhancement in the production of high-£'T pho- 
tons [l3|- We take advantage of this enhancement by 
comparing the photon Et distribution in data to both 
SM and beyond-the-SM predictions. Binned maximum 
likelihood measurements of the coupling parameters that 
describe Z7 interactions in the Lagrangian are performed. 



We calculate separate likelihoods for the Z l^l^ and 
Z ^ vv samples and combine the likelihoods to produce 
the final result. 

The CDF detector is covered in detail elsewhere [3, 
15| . The transverse momenta (j)t) of charged parti- 
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cles are measured by an eight-layer silicon strip detec- 
tor [3l and a 96-laycr drift chamber (COT) 17 1 inside 
a 1.4 T magnetic field. The COT provides tracking cov- 
erage with high efficiency for the pseudorapidity range 
I77I < 1 [l3|- Electromagnetic and hadronic calorimeters 
surround the tracking system. They arc segmented in a 
projective tower geometry and measure the energies of 
charged and neutral particles in the central (|?y| < 1.1) 
and forward (1.1 < \r\\ < 3.6) regions. Each calorimeter 
has an electromagnetic shower profile detector positioned 
at the shower maximum [l8l |. The calorimeters arc sur- 
rounded by drift chambers that detect muons. 

The measurements of anomalous trilinear gauge cou- 
pling parameters in the Z — >■ and Z ^ vv decay 
channels differ both in event selection and background 
estimation. For the Z7 — > decay channel we iden- 
tify events containing Z — )■ fJ-^ fJ.~ and Z — e~^e~ candi- 
dates along with prompt photon candidates with E^ > 50 
GeV. According to experiments performed on simulated 
events, this choice of E^ requirement maximizes the abil- 
ity of the analysis to exclude anomalous couplings as- 
suming SM physics, although a serious loss in sensitiv- 
ity only occurs if the E^ requirement is placed at 100 
GeV or higher. The previous CDF analysis used a much 
less restrictive requirement of E^ > 7 GeV, as the Z7 
cross section was being measured in addition to trilinear 
gauge coupling parameters Q; additionally, placing the 
cut at 50 GeV allows for a control region to be based off 
of lower- photons. Event selection starts with inclu- 
sive muon (electron) triggers that require muon > 18 
GeV/c (electron Et > 18 GeV). For electrons, a track 
must be reconstructed in the COT or in the silicon de- 
tector; additionally, the energy deposited by the candi- 
date in the calorimeter must be isolated. For muons, a 
track must be reconstructed in the COT; additionally, no 
more than a few GeV of energy may be deposited in the 
calorimeters so that the candidate is compatible with a 
minimum ionizing particle. The two lepton candidates 
must correspond to the same flavor, with a requirement 
of pt > 20 GeV/c (Et > 20 GeV) on one muon (elec- 
tron) candidate and pT > 10 GeV/c {Et > 10 GcV) on 
the other; furthermore, if the charges of both leptons are 
well-measured, the signs of these charges must be oppo- 
site. Studies of the invariant mass distributions of the 
two lepton candidates indicate that we retain a very high 
purity of Z bosons (over 99%) despite the loose selection 
requirements. 

Once we have selected events with Z l^l^ can- 
didates, we look for isolated photons that pass stan- 
dard CDF requirements [l9| in the central region (I77I < 
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1.1) with > 50 GcV and are well separated 

from the Z decay leptons (Ai?^^ > 0.7, with AR = 

^ — (jy-yY + {vi ~ ViY)- Additionally, we require that 
the two lepton candidates and the photon candidate form 
a three-body invariant mass greater than 100 GeV/c^ in 
order to discriminate against events where the photon is 
radiated from one of the leptons from the Z boson decay. 
The estimated contribution of SM Z^ events is derived 
from Monte Carlo (MC) simulations that use the Baur- 
Berger package at the generator level [l3| and PYTHIA [2^ 
for particle showering. This method yields a prediction of 
87.2 ± 7.8 Z^ events that pass our selection requirements, 
where the uncertainty is dominated by the uncertainty on 
the luminosity and the predicted cross section. The non- 
Z^ events that pass these selection requirements result 
from hadronic jets being reconstructed as prompt pho- 
tons and leptons (more commonly electrons). This back- 
ground is estimated by calculating separate probabilities 
for a jet to mimic a photon or lepton as a function of 
jet Et, and applying them to jets in events to which all 
our requirements have been applied except those pertain- 
ing to the mimicked particle. For photons and electrons, 
these probabilities are calculated by taking the ratio of 
the number of individual photon or electron candidates 
to the number of jets in a sample of data events where 
only the presence of at least one jet is required. The num- 
ber of photon and electron candidates is corrected for the 
expected contribution of true photons or electrons in this 
sample. We estimate the probability for a false muon can- 
didate from the number of dimuon Z decay candidates in 
which both muon candidates have the same charge. Over- 
all, the non-Zj background contribution is very low: of 
the 91 events that pass our requirements, less than one 
event involving a mimicked photon or lepton is expected. 

In order to identify Zj candidate events in the Z ^ vv 
decay channel, we require solitary high- i?T photons and 
a transverse energy imbalance [21| in the detector. These 
events must pass a trigger requirement of an electromag- 
netic cluster with Et > 25 GeV and jryj < 1.1 as well as 
missing transverse energy in excess of 25 GeV. For our 
signal region we require Ej. > 100 GeV, a threshold op- 
timized in the same manner as the Z — ^ case. To 
account for the neutrinos we require a transverse energy 
imbalance of at least 50 GeV. In order to discriminate 
against W boson contamination in our sample, we re- 
ject events containing any tracks with px > 10 GeV, any 
electron candidates with Et > 15 GeV, or any muon 
candidates with pT > 10 GcV/c. Additionally, we reject 
events that have any jets with Et > 15 GeV in order 
to reduce the mismeasurcmcnt of missing transverse en- 
ergy. The primary SM source for photons passing these 
requirements is Z7 events in which the Z has decayed to 
a pair of neutrinos, as shown in Table ID The method of 
estimating the expected number of Zj events is the same 
as that used for the Z — >■ candidate sample. 



The primary source of non-Zj events in the final 
Z"f — > i^Dj candidate sample is cosmic ray interactions. 
High-Ex photons from cosmic rays leave large transverse 
energy imbalances in our detector, mimicking the pres- 
ence of neutrinos. Therefore, additional event require- 
ments are applied to reduce the contributions from cos- 
mic ray events. First, we require that the energy de- 
posited in the electromagnetic calorimeter appear within 
a timing window centered on the pp interaction. Second, 
we use a relevance vector machine (RVM) multivariate 
discriminator [22| to distinguish whether a photon came 
from a collision or a non-collision source; the three inputs 
used for the RVM discriminator are the 0-angle between 
the photon candidate and the closest muon candidate (if 
any) , the ratio of energies from the photon candidate in 
the electromagnetic and hadronic calorimeters, and the 
ratio of energies from the electromagnetic shower profile 
detector and the electromagnetic calorimeter. We use 
photons outside the timing window to train the RVM for 
non-collision sources, and photons recoiling against jets 
to train for collision sources. The RVM discriminator 
reduces the contribution from cosmic ray events by an 
additional 90%. Finally, we require the event to have a 
reconstructed vertex of at least three tracks from a pp 
interaction. After applying these selection requirements, 
we have 85 candidate events in our sample. Despite the 
anti-cosmic ray requirements, cosmic ray events remain 
the second largest contributor to our sample, after Z7 
events. 

We model two other major categories of non-Zj events: 
one in which a charged lepton from W — ^ ev, W — >■ fiv, 
or W ^ TV decay is reconstructed as a photon, and the 
other in which a true photon is produced but another ob- 
ject (e.g. a lepton) is lost or only partially reconstructed, 
creating a large transverse energy imbalance. For the for- 
mer case, the rate at which electrons are reconstructed as 
photons in the detector has been calculated using events 
with an electron and photon pair candidate that has an 
invariant mass near the mass of the Z, i.e., events in 
which the photon candidates are almost entirely elec- 
trons in actuality. The rate at which ^'s and r's are 
reconstructed as photons is taken from MC. For the lat- 
ter case, which encompasses — > Ivj events in which 
a lepton is lost and 77 events in which a photon is lost, a 
two-step process is used to calculate the expected number 
of events. First, events in data are selected such that we 
obtain a very pure sample of one of the aforementioned 
event types in which there is no lost object. Then, we cal- 
culate the fraction of the corresponding events in MC in 
which an object is not detected, and this fraction is used 
to scale the photon Et distribution of the data events 
so as to provide an estimate of this background's photon 
Et distribution in the signal sample. An exception to 
this method is the case in which a r is lost; due to the 
difficulty of reliably identifying r candidates, this back- 
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ground is estimated purely from MC simulations. Fur- 
ther details on these methods of background prediction 
can be found in [2^, a CDF analysis which used very 
similar event requirements. We see excellent agreement 
between the SM predictions and the data in the control 
regions of 15 < < 40 GcV {Z case) and 

70 < < 100 GeV [Z -> vv case). 



Process 


Events 


Z7 — > vv^ 


52.8 ±4.6 


cosmics 


14.9 ± 1.4 


W ^ ev 


3.9 ±0.8 


W — ^/tv^ 


1.6 ±0.3 


Wj — ?> 61^7 


1.1± 1.1 


W"f — > fii>^ 


1.8± 1.3 




4.5 ± 1.3 


77 


5.3 ± 1.9 


SM Total 


85.9 ±5.6 


Data 


85 



TABLE I: SM expected contributions to the Z'j vv'y can- 
didate sample. Uncertainties shown are systematic only and 
thus exclude the statistical uncertainties. 



Assuming gauge and Lorentz invariance, eight parame- 
ters are needed to describe couplings, denoted by hj^ 
where V is either a or a 7 and the index i runs from 1 
to 4; these parameters are all zero at tree level [T^]. Inter- 
action amplitudes are linear in these parameters. Indices 
1 and 2 represent CP-violating terms while indices 3 and 
4 represent CP-conserving terms. We assume CP conser- 
vation in these interactions by setting hX^ = h^^ = and 
we investigate the possibility of non-zero values for /i^ 
and /i4^, corresponding to electric dipole and magnetic 
quadrupole transition moments 2J] . In order to preserve 
unitarity at large incoming parton center-of-mass energy 
an s-dependent form factor is used to suppress the 
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coupling, constructed as h^{s) ~ (i+s/A^)" ' ■^'^ci'c n — i 
for hX^Q and /140 [l^l ■ The parameter A describes the pre- 
dicted energy scale of the new physics that creates anoma- 
lous Z7 couplings. 

For a given set of anomalous coupling parameter values, 
we compute a likelihood for the i?^ distribution. Hence, 
we have Jlj^^i -^(^il^'-D' where Xj represents the number 
of entries in the jth of N bins in our Elp distribution 
and }i( denotes the coupling parameter being measured 
(the other three being held fixed at zero). The bin- by- 
bin likelihood L is simply the Poisson probability of the 
number of observed entries given the expected number 
of entries for the value of }i( . This limit method re- 
quires a predicted Et distribution for each combination 
of the four coupling parameters. To create these dis- 
tributions, we produce Z7 MC events at the generator 
level using the Baur-Berger package [l3|- Modeling the 



particle showering process and detector response in MC 
separately for every parameter value is computationally 
impractical. To mimic fully-simulated MC events we first 
determine the efficiency for a generated event to pass all 
of the event requirements as a function of generator-level 
i?J and \rp\; these functions are derived from a SM MC 
sample which has used the full simulation of the detector. 
Due to the correlation between El^ and the Z kinematics, 
we create and combine separate templates for the cases 
of central-central, central-forward, and forward-forward 
lepton pairs, "central" denoting < \t']\ < 1.1 and "for- 
ward" denoting 1.1 < \r]\ < 2.8. We then apply this 
efficiency function to generator-level MC samples to get 
the expected E^ distributions. The final prediction is 
the sum of this Zj prediction with the predictions of the 
non-Zj backgrounds. 

In Fig.lll for both the Z l^l^ and Z ^ vv cases, the 
Elp distributions in data are compared to the SM predic- 
tion and beyond-tlic-SM predictions; it can be seen that 
the production of high-£^T photons is far more likely in 
the beyond-the-SM cases compared to the SM case. The 
uncertainty bands shown for the SM predictions illustrate 
the systematic uncertainties on those predictions. These 
uncertainties are dominated by the 7% uncertainty on 
the theoretical Z7 cross section [251 and the 6% uncer- 
tainty on the luminosity [26]; the other sources are the 
reconstructed photon's energy scale and efficiency, as well 
as uncertainties on the number of non-Zj background 
events. The effect of these systematic uncertainties on 
the limits is negligible — of the order of a couple of per- 
cent of the limit values. 

With the likelihood distribution for a given hY , taking 
a flat Bayesian prior in hY allows us to set Bayesian credi- 
bility limits on the parameter. These limits are defined as 
the values of hY which demarcate the central 95% of the 
integral of the likelihood distribution. The resulting al- 
lowed ranges for the strength of anomalous couplings are 
shown in Table HI The values A = 1.2 TeV and A = 1.5 
TeV have been chosen to allow direct comparisons with 
earlier CDF [l[ and DO 0| results, respectively. We see 
no evidence for anomalous couplings. 

In conclusion, 

tons produced in association with Z bosons in both the 
Z ^ vv and Z — > decay channels in a data sam- 

ple corresponding to an integrated luminosity of approxi- 
mately 5 fb~^ is consistent with SM couplings. We place 
95% Bayesian credibility limits of \hCl'^\ < 0.027 and 
\hj'^\ < 0.0013 on the CP-conserving Z7 couplings at 
A = 1.2 TeV and \hl^\ < 0.022 and \hl'^\ < 0.0009 at 
A = 1.5 TeV; these are significantly tighter constraints 
on beyond-thc-SM contributions than those provided by 
previously measured limits. 
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Zy^tl'y Data 

Standard Model Prediction 
= 0.032 Prediction 
= 0.0013 Prediction 



Photon (GeV) 



10=^ 



Parameter 


(A = 1.2 TeV) 


(A = 1.5 TeV) 




-0.024, 0.027 


-0.020, 0.021 


hi 


-0.0013, 0.0013 


-0.0009, 0.0009 


hi 


-0.026, 0.026 


-0.022, 0.020 


hi 


-0.0012, 0.0013 


-0.0008, 0.0008 



TABLE II: Allowed ranges (95% Bayesian credibility limits) of 
anomalous couplings for A = 1.2 and 1.5 TeV using nota- 
tion from Ref. jT^]. Each parameter's limits are set assuming 
the other three parameters have values fixed at 0. 
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FIG. 1: Comparison of the measured ET^ distribution with the 
predicted distributions from both the SM and beyond-the-SM 
scenarios for Z — >■ V^l' (top) and Z ^ vv (bottom) candi- 
date samples, at A = 1.5 TeV. The beyond-the-SM scenarios 
chosen here can be excluded at 95% Bayesian credibility level 
in each sample. Note the greatest difference in the SM and 
beyond-the-SM cases is found offscale at E~Jp> 200 GeV; the 
lack of data events in this region indicates good agreement 
with the SM. 
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